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Following the early prediction of the skyrmion 
lattice (SkL)33®—a periodic array of spin 
vortices—it has been observed recently in vari¬ 
ous magnetic crystals mostly with chiral struc¬ 
ture^!^! Although non-chiral but polar crys¬ 
tals with G nv symmetry were identified as ideal 
SkL hosts in pioneering theoretical studies 1 4 , 
this archetype of SkL has remained experimen¬ 
tally unexplored. Here, we report the discov¬ 
ery of a SkL in the polar magnetic semiconductor 
GaV^Sg with rhombohedral (C 3 V ) symmetry and 
easy axis anisotropy. The SkL exists over an un¬ 
usually broad temperature range compared with 
other bulk crystals 6 12 and the orientation of the 
vortices is not controlled by the external mag¬ 
netic field but instead confined to the magnetic 
easy axis. Supporting theory attributes these 
unique features to a new non-chiral or Neel-type 
of SkL 1 2 24 describable as a superposition of spin 
cycloids in contrast to the Bloch-type SkL in chi¬ 
ral magnets described in terms of spin helices 6 17 . 

In non-centrosymmetric crystals, the energy associated 
with ferromagnetic domain walls can be negative, hence, 
the homogeneous ferro magn etic state becomes unstable 
against SkL formation 1 5 25 . Two basic types of mag¬ 
netic domain walls can form different skyrmionic spin 
textures 1 3 2 3 Bloch-type domain walls, where the spins 
rotate in the plane parallel to the domain boundary, can 
form whirlpool-like skyrmions with a given handedness 
as observed in chiral magnets 6 17 . In contrast, Neel-type 
domain walls, where the spins rotate in a plane per¬ 
pendicular to the domain boundary, can produce non- 
chiral skyrmions. (See Fig. 5 for an illustration of the 
two skyrmion types.) This second skyrmion archetype, 
expected t o em erge in polar magnets with C nv crystal 
symmetr}®^!], has not been observed yet. We reveal 
the formation of such a SkL in GaV 4 S 8 , a magnetic semi¬ 
conductor from this crystal symmetry class, via magnetic 


susceptibility, atomic force microscopy (AFM) and small- 
angle neutron scattering (SANS) measurements. 

GaA^Sg, a member of the lacunar spinel famil})^^, 
has a non-centrosymmetric cubic (T^) structure at room 
temperature 32 . Lacunar describes the lack of every sec¬ 
ond Ga atom compared to the normal spinel structure. 
These ordered defects break the vanadium pyrochlore 
lattice into a network of alternating larger and smaller 
V 4 tetrahedra, the so-called breathing pyrochlore lattice. 
The magnetic building blocks are the smaller V 4 clusters 
with spin 1 / 2 , and form a face centered cubic (FCC) lat¬ 
tice. The hybridization on a single V 4 unit leads to one 
unpaired electron occupying a triply degenerate cluster 
orbital. This orbital degeneracy is lifted by a coopera¬ 
tive Jahn-Teller distortion that drives a cubic to rhombo¬ 
hedral structural transition at T S =42K 32 . In the polar 
rhombohedral (C 3 V ) phase the FCC lattice is stretched 
along one of the four cubic ( 111 ) axes 32 . 

The structural transition creates a multi-domain state 
with submicron-thick sheets of the four different rhom¬ 
bohedral domains, as seen in Fig. la. At T s the magnetic 
exchange interaction changes from antiferro- to ferromag¬ 
netic and the material undergoes a magnetic transition 
at Tc —13 K, a temperature very close to the Curie-Weiss 
temperature in the rhombohedral phasd^. In the inset of 
Fig. lb, the dependence of the magnetization curves on 
the field orientation indicate a considerable anisotropy. 
Indeed, our mean-field analysis shows the ground state to 
be an easy axis ferromagnet with Jy^8.7 K and J±&8.3 K 
as described in the Supplementary Information. Here, 
J|| and J± respectively denote the exchange coupling for 
spin components parallel and perpendicular to the easy 
axis, which itself coincides with the direction of the rhom¬ 
bohedral distortion. These exchange values are also con¬ 
sistent with Tc =13 K. 

GNV^Sg was considered as an ordinary ferromagnet 
below Xd 32 33 l though there is a report about low-field 
magnetic anomalies in this material-^. We found that 
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FIG. 1: | Magnetic phases in the lacunar spinel GaV^Sg. a, Topographic image recorded by AFM on the (100) surface 
of the crystal at T— 8.9K. The colour scale corresponds to the altitude perpendicular to the image plane. Alternating blue 
and yellow stripes running along the cubic [001] direction are structural domains with different rhombohedral axes, while edges 
along the [010] axis are epilayer terraces. The yellow and white dashed lines highlight a domain boundary and a step between 
epilayer terraces, respectively. Magnetic patterns in Fig. 2 were recorded within the area of the red square, b, Magnetization 
curves measured at 12 K in B || [100], [110] and [111] (shifted vertically for clarity). Magnetization steps are observed at different 
field values dependent on the orientation. The inset shows the magnetization measured at T— 2 K up to B— 1.2 T. d-f, Magnetic 
phase diagrams derived from the field- and the temperature-dependence of the magnetization. Circles and squares correspond 
to peaks in the field- and temperature-derivative of the magnetization curves, respectively. The insets show the orientation 
of the magnetic field relative to the easy axes of the four rhombohedral domains (cubic body diagonals). The easy axes of 
magnetically favoured/unfavoured domains are indicated by thick/thin green lines. For B||[111] , besides the cycloidal and the 
SkL states, there are two additional phases extending up to higher fields labelled as the cycloidal* and the SkL* states, c, 
Phase boundaries from panels d-f, after projecting the magnetic field onto the easy axis of the corresponding domains. 


the magnetic state is far from ordinary and contains 
several neighbouring phases. Below Tc, the low-field 
magnetization displays a sequence of steps indicating 
metamagnetic transitions (see Fig. lb). The positions 
of the steps, associated with sharp peaks in the field- 
derivative of the magnetization, are different for mag¬ 
netic field, B, applied parallel to the [111], [110] and 
[100] cubic axes. The corresponding phase diagrams 
are shown in Figs. ld-f. Similar low-field steps have 
been observed to sep arate the helical, SkL and conical 
states in MnSP™, Mn 1 _ a .(Co,Fe) a .SP3, FeG<P and 
Cu 2 OSeOP. In GaV 4 S 8 we assign the three differ¬ 
ent magnetic phases as cycloidal, SkL and ferromagnetic 
based on AFM imaging, SANS experiments and theoret¬ 
ical calculations as described below. 

For B|| [111] there exist two additional magnetic states 


besides the cycloidal, SkL and ferromagnetic states. This 
arises due to multi-domain nature of the crystal and the 
easy axis anisotropy, since only one of the domains has its 
easy axis parallel to B, while the easy axes of the other 
three domains—[111], [111] and [Ill]—lie at 71 degrees 
to B. For these latter three domains the non-collinear 
spin structures survive to higher fields because the mag¬ 
netic anisotropy energy is much larger than the Zeeman 
energy, and only the field component parallel to the easy 
axis of each domain can influence its magnetic state. 
This leads to the two additional phases persisting up to 
5=80 mT and 160 mT, marked with asterisks in Fig. Id, 
and assigned as the cycloidal and SkL states common for 
these three domains. For B||[110], the cycloidal and SkL 
phases common for the two domains with easy axes [111] 
and [111] are indicated in Fig. le. The other two domains 
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FIG. 2: | Real-space imaging of the magnetic patterns in GaXGSs. a, AFM image recorded at T— 8.5K and in 
5=20 mT on the (100) surface in the area indicated by the red square in Fig. la. Colour coding corresponds to dissipated 
power due to magnetic interactions between the tip and the sample. A modulated structure with a single q- vector and a 
corresponding periodicity of a cyc =17.7±0.4 nm is observed on the left side of the domain boundary (yellow dashed line). A 
similar pattern rotated by about 90 degrees emerges on the other domain, though the magnetic contrast is weaker, b, Magnetic 
pattern measured in the same area at T— 8.9K and in 5=50 mT. Distorted triangular lattices, rotated relative to each other 
by approximately 90 degrees, are observed in both domains with better contrast on the right side. The area with red border 
is scaled by a/ 3 according to the red arrows, i.e. perpendicular to the direction of the shortest lattice periodicity. The area, 
transformed to a square shape after this contraction, is displayed in panel c. The green spots indicate a regular triangular SkL 
with a lattice constant of a s k y =22.2±l nm and vortex cores parallel to the [111] axis. For more details about the data analysis 
and the evolution of the magnetic pattern as a function of magnetic field see the Supplementary Information, d, Two sections 
of a hexagonal array of vortex lines running along the [111] axis. The intersections of the vortices with the (111) and (100) 
planes respectively form a regular and distorted triangular lattice. 


with easy axes perpendicular to B cannot contribute to 
the magnetization is this field range. For B|| [100], the 
easy axis of each domain span the same angle to B, hence 
all domains share common phase boundaries as shown in 
Fig. If. Phase boundaries for all field orientations can be 
scaled together by projecting the field onto the easy axes 
of the corresponding domains (see Fig. lc). Surprisingly, 
the SkL state exists over an unusually broad tempera¬ 
ture range, down to ^0.68Xc, compared with other bulk 
crystals where this phase is stable only over a 2-3% range 
immediately below 

To observe directly the spin pattern in real space we 
performed AFM imaging in magnetic fields applied per¬ 
pendicular to the sample surface using magnetic can¬ 
tilever tipsP. For B ||[100], Figs. 2a and b show images 
of the same surface area taken at T=8.5K and 8.9 K in 
5=20 mT and 50 mT, respectively. The trail running 
along the [001] axis and dividing each image into two 
halves is a structural domain boundary. 

Although all four domains should host the same mag¬ 
netic states for B|| [100], we found that the magnetic pat¬ 
tern always changed at the domain boundaries, and often 
it was not possible to achieve good contrast for the two 
sides simultaneously. Since the wave vectors of the mag¬ 
netic modulation (q-vector) are expected to be perpen¬ 
dicular to the easy axis, they can be different for different 
domains, resulting in a rotation of the magnetic pattern 
between the two sides. In Fig. 2a, a modulated structure 
with a single q-vector, representative of the zero- and 


low-field region, is observed with good contrast on the 
left side. Our calculations described later reproduce this 
phase as a cycloidal spin state, where the spins rotate in 
a plane containing the q-vector. This is in contrast to 
the helical order generally found in other SkL materials, 
where the spins rotate in the plane perpendicular to q. 
Assuming that the q-vector is perpendicular to the easy 
axis, which is either the [111] or [111] axis for the left¬ 
side domain, the periodicity of the cycloid in 5=20 mT 
is a cyc =17.7=t0.4nm. (For more details see the Supple¬ 
mentary Information.) This value is confirmed by our 
SANS study. 

With increasing field the cycloidal pattern transforms 
into a distorted triangular lattice as shown in Fig. 2b. 
The same pattern rotated by about 90 degrees appears 
on both sides of the domain wall, though the contrast 
is weaker on the left side. As demonstrated in Fig. 2c, 
re-scaling the lattice on the right side by \/3 perpen¬ 
dicular to the direction of the shortest lattice constant 
leads to a regular triangular lattice with a lattice con¬ 
stant of a s k y =22.2±l nm. This value is also consistent 
with our SANS data. The observed distortion is nat¬ 
urally explained by the orientational confinement of the 
SkL. Since the vortex cores are not aligned parallel to the 
magnetic field but instead the easy axis, which is either 
the [111] or [Ill] axis for the right-side domain, they form 
a regular triangular lattice in the hard plane. Therefore, 
the terminus of the vortices on the (100) surface form a 
distorted triangular lattice as sketched in Fig. 2d. When 
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FIG. 3: | Small angle neutron scattering study of the magnetic states in GaV^Ss. a-d, SANS images measured 
at T=10.9K in various fields for B||[111]. The images were recorded at field values representative of the different magnetic 
states assigned in Fig. Id. e &; f , SANS images taken at T=10.9K in different magnetic fields for B|| [100]. Labels refer to the 
magnetic states assigned in Fig. If. The dashed circles having common diameter in all the images help to visualize the change 
in the magnitude of the q-vectors. The two colour bars indicate the scattering intensity for the two field orientations, g &; 
h, The magnetic periodicity, a = 27r/|q|, (left scales) and the total scattering intensity (right scales) as a function of the field 
strength for the two magnetic field orientations. Vertical dashed lines indicate the magnetic phase boundaries according to 
Fig. Id and f. The splitting of each Bragg peak into two peaks with different |q| was resolved for B|| [111] in the field range of 
F>=44-64mT (see the Supplementary Information) which arises from the coexistence of cycloidal* and SkL states in different 
domains. 


further increasing the field, the magnetic pattern disap¬ 
pears indicating a uniform ferromagnetic state. 

SANS was employed to confirm the bulk nature of the 
magnetic structures and determine the q-vectors directly. 
Sharp magnetic Bragg peaks were always observed at q- 
vectors parallel to the cubic (110) directions irrespective 
of the field orientation, as shown in Fig. 3. For B||[111] 
the ±q-vectors form a regular hexagon, while for B|| [100] 
a square of four spots is observed. These results, in agree¬ 
ment with the results of AFM imaging, confirm the ori¬ 
entational confinement of the SkL, i.e. the alignment of 
the vortex lines parallel to the easy axis of each domain. 
It also shows that in both the cycloidal and SkL states, 
the three possible ±q pairs are parallel to the three [110] 
axes lying in the hard plane. For B||[111] a pale ring 
structure appears besides the hexagon in low fields. This 
implies that the q-vectors are not fully pinned to the [110] 
axes but that they can point in any direction within the 
hard plane, at least in the surface and domain boundary 
regions observed by AFM. 

In each rhombohedral domain the cycloidal state can 
emerge with six different q-vectors, while the SkL state 
is unique, and describable as a superposition of three cy¬ 
cloids whose q-vectors sum to zero. First-order Bragg 
peaks observed using SANS do not allow a direct dis¬ 


tinction between the cycloidal and the SkL states due to 
the presence of magnetic and structural domains. Never¬ 
theless, the field dependence of the magnetic periodicity 
(a = 27r/|q|) clearly indicates the phase boundaries for 
both field orientations as shown in Fig. 3h. The coex¬ 
istence of the cycloidal* and SkL states in different do¬ 
mains for B||[111] is also traced via the splitting of |q|. 

In GaV 4 S 8 , the S=l/2 spins sit on an FCC lattice 
stretched along the [111] direction, which can be regarded 
as triangular lattices stacked along [111]. To determine 
the spin patterns in the different phases we studied the 
following classical Heisenberg model on the triangular 
lattice by a Monte Carlo technique; 

^ = ^ex T HdM T ^Zeeman? (i) 

where 

T^ex — ^ ^ A_WlxiWlxj T JA_WlyiWlyj T ? 

<i,j> 

"Hdm = ^2 D ij( m i xrrij), 

<i,j> 

^Zeeman — B \ ^ Tfl z i • 
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FIG. 4: | Spin patterns in the magnetic phases of GaV^Ss. a, FCC lattice of V 4 units each carrying a spin 1/2 and the 
orientation of the Dzyaloshinskii-Moriya vectors for bonds on the triangular lattice within the (111) plane (chosen as the xy 
plane in the calculation), b, Cycloidal spin state obtained for the spin model in Eq. 0 on the triangular lattice, c, Magnified 
view of the magnetization configuration for the cycloidal state, d, Bragg peaks (q-vectors) of the cycloidal state in panel b in 
reciprocal space, e, SkL state obtained for the spin model in Eq. ([!]) on the triangular lattice, f, The magnified view of the 
magnetization configuration for the SkL state clearly shows the Neel-type domain wall alignment, g, Bragg peaks of the SkL 
state in panel e. The q-vectors of the Bragg peaks are located along the [110] directions (white lines) in the hard plane, for 
both the cycloidal and SkL states. 


For the exchange interaction 7i ex , we consider ferromag¬ 
netic coupling with XXZ-type anisotropy, where the z 
(easy) axis is parallel to the [111] axis. The direction of 
the Dzyaloshinskii-Moriya vectors, Dij , are determined 
by the FCC crystal symmetry as shown in Fig. [4 ^l. For 
the calculation we chose J||/J_l= 1.06 and D/J_l= 0.35, 
where D = | D ^ |. 

This model provides three ordered phases: the cy¬ 
cloidal, SkL and ferromagnetic states. Figure 4 dis¬ 
plays the spin patterns obtained for the cycloidal and 
SkL states. The skyrmions exhibit the Neel-type domain 
wall alignment along the radial direction from their cores 
to their peripheries. Namely, the spins rotate within the 
plane parallel to the radial direction. The Fourier compo¬ 
nents of the magnetization configurations in the cycloidal 
and the SkL states are displayed in Figs. 4d and g, re¬ 
spectively. These results, in full agreement with SANS 
data, show that the ±q- vectors are parallel to the [110] 
axes in the hard plane in both states. 

Extension of the model in Eq. 0 to the distorted FCC 
lattice requires additional parameters as both J\\/J± and 
D/J± can be different for bonds within and between the 
triangular lattice planes. Instead, we adapted the contin¬ 
uum model developed by Bogdanov and Hubert for C nv 
symmetr^ to the cubic lattice in order to reproduce the 
main features of the phase diagram. Here, s are per¬ 
pendicular to the bonds along the x and y directions as 
shown in the inset of Fig. 5c, and values of J||/J_l=1.06 
and D/J_l=0.48 were used. We stress that irrespective 
of the details of possible models, GaV^Ss can only host 
a Neel-type SkL due to the orientation of Dzyaloshinskii- 


Moriya vectors governed by its C% v symmetrjB®^!. 

Figure 5c displays the phase diagram obtained for the 
simple cubic lattice. In agreement with experiment the 
triangular SkL appears only at finite temperatures. The 
skyrmions have the same Neel-type pattern, shown in 
Fig. 5b, as obtained for the triangular lattice model. 
The cycloidal phase remains stable even at T=0 in con¬ 
trast to the experimental observation. This discrepancy 
likely arises from not including corrections due to quan¬ 
tum fluctuations in our classical Heisenberg model. 

In a new class of materials we observed a SkL 
state emerging over an extraordinarily broad tempera¬ 
ture range. While chiral or Bloch-type skyrmions have 
been investigated extensively in chiral magnets, lacunar 
spinels provide a unique arena to study another topolog¬ 
ical spin pattern, the non-chiral or Neel-type skyrmions. 
We found that the strong orientational confinement of the 
vortices ensures the robustness of two distinct skyrmionic 
states with a core magnetization pointing either up or 
down the easy axis. This may facilitate a unique mag¬ 
netic control of the SkL, since we expect that the SkL in¬ 
duced by a magnetic field parallel to the easy axis can be 
rotated within the hard plane by an additional transverse 
field component. Such magnetic control is inconceivable 
in cubic helimagnets where vortex cores instantaneously 
co-align with the magnetic field. In addition, the polar 
crystal structure of lacunar spinels 31 32 may be exploited 
for a non-dissipative electric field control of the SkL. 
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FIG. 5: | Neel-type SkL in GaXGSs. a, The Bloch-type 
skyrmion is a chiral object, i.e. it has no mirror-plane symme¬ 
try. b, In contrast, the Neel-type skyrmion has mirror-plane 
symmetry (combined with time reversal) for vertical planes 
through its center, c, Theoretically obtained temperature (T) 
versus magnetic field (B) phase diagram for the spin model 
by Eq. Q on the cubic lattice. The inset shows the direction 
of the Dzyaloshinskii-Moriya vectors for bonds within the xy 
plane. The skyrmion lattice state corresponds to a Neel-type 
SkL. 


Methods 


Sample synthesis and characterization. Single 
crystals of GaA^Ss were grown by the chemical vapour 
transport method using iodine as the transport agent. 
The crystals are typically cuboids or hexagonal slabs with 
masses ranging from 1-60 mg. The sample quality was 
checked by powder X-ray diffraction, specific heat and 
magnetization measurements. The crystallographic ori¬ 
entation of the samples was determined by X-ray Laue 
and/or neutron diffraction prior to the magnetization, 
SANS and AFM studies. The magnetization measure¬ 
ments were performed using an MPMS from Quantum 
Design. The field dependence of the magnetization was 
measured in increasing temperature steps following an 
initial zero-field cooling to 2K. The temperature depen¬ 
dence of the magnetization was also measured in different 
fields. For both of these approaches, we found no hys¬ 
teresis in the magnetization curves thus indicating that 
the SkL state cannot be quenched, i.e. stabilized as a 
metastable state, in this compound. 

Non-contact atomic force microscopy. All AFM 
data presented were obtained by non-contact atomic 
force microscopy (nc-AFM) performed in an Omicron 


cryogenic ultra-high vacuum (UHV) STM/AFM instru¬ 
ment using the RHK R9-control electronics. The micro¬ 
scope features interferometric detection which as a side- 
effect allows the tuning of the cantilevers’ effective qual¬ 
ity factor Q e ff by varying the laser intensity used for 
the detection of the cantilever motion. For all measure¬ 
ments we used magnetically coated AFM-tips from the 
SSS-QMFMR series from NANOSENSORS™in order to 
achieve sensitivity to magnetic forces. 

In nc-AFM mode the frequency / of an oscillating can¬ 
tilever is measured via a phase locked loop and kept con¬ 
stant by the topography feedback loop. The frequency 
shift A/ with respect to the resonance frequency /o far 
away from the sample provides a measure of the force 
gradient apparent between tip and sample. Additionally, 
the oscillation amplitude A is kept constant by adjusting 
the excitation amplitude A exc . Electrostatic interactions 
between tip and sample have been minimized using an ad¬ 
ditional FM-Kelvin control loop without any additional 
lateral information in the measured contact potential dif¬ 
ference. 

We could not obtain magnetic contrast in magnetic 
force microscopy mode, where the tip is retracted from 
the surface and the conservative magnetic forces are mea¬ 
sured via the induced frequency shift A/. Instead, a 
dissipative force interaction induced the magnetic con¬ 
trast in the excitation channel when the tip was brought 
close to the sample surface. The nature of this dissipa¬ 
tive interaction is subject to further investigations. We 
could obtain sufficient magnetic contrast only in images 
recorded below 10 K, which is likely related to the low 
spin density of the material corresponding to a single spin 
per formula unit and the strong fluctuations of the S=l/2 
spins. 

With knowledge of the cantilever’s spring constant &, 
the actual oscillation amplitude A , the resonance fre¬ 
quency /, the effective quality factor Q e ff,o and the 
signal-to-drive ratio Rq far away from the sample we com¬ 
pute the dissipated power from the measured excitation 
amplitude as 


irkfRoA 

Qeff, 


For the images in figure 2 these parameters are k = 
3.4 N/m, A = 3nm, Q e//?0 = 4.28 • 10 5 , / = 70.87kHz, 
Ro = 30.92 /im/V with excitation amplitudes ranging 
from 150/W to 1.1 mV. 

Small-angle neutron scattering. Small-angle neu¬ 
tron scattering (SANS) was used to study the long- 
wavelength, microscopic magnetic states in a 25.5 mg sin¬ 
gle crystal sample of GaV 4 S 8 - The SANS measurements 
were carried out using the D33 instrument at the Institut 
Laue-Langevin (ILL), Grenoble, France, and the SANS- 
II instrument at the Swiss Spallation Neutron Source 
(SINQ), Paul Scherrer Institut, Switzerland. In a typi¬ 
cal instrument configuration, neutrons of wavelength 5 A 
were selected with a FWHM spread (AA/A) of 10%, 
and collimated over a distance of 5.3 m before the sam- 
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pie. The scattered neutrons were collected by a two- 
dimensional multi-detector placed 5 m behind the sam¬ 
ple. 

The single crystal sample of GaA^Ss was mounted in¬ 
side a horizontal field cryomagnet that was installed on 
the beamline so that the direction of applied magnetic 
field was approximately parallel to the neutron beam. 
The sample was oriented to have a horizontal scatter¬ 
ing plane that included the [100] and [111] crystal direc¬ 
tions. Applying the magnetic field along either of these 
directions was achieved by simply rotating the sample 
stick inside the cryomagnet. The SANS measurements 
were done by both tilting and rotating the sample and 
cryomagnet together through angular ranges that moved 
the magnetic diffraction peaks through the Ewald sphere. 
The SANS patterns presented in Fig. 3 were constructed 
by summing together the detector measurements taken 
at each rotation angle, which allows all of the Bragg spots 
to be presented in a single image. 

Theoretical calculations. The Monte-Carlo calcu¬ 
lations for Fig. [4] are performed for the model in Eq. [T| 
on the two-dimensional triangular lattice with 144 x 144 
sites, while those for Fig. [5] are for the model in Eq. [T| 
on the three-dimensional cubic lattice with 72 x 72 x 8 
sites. For both cases periodic boundary conditions are 
imposed. The XXZ-type anisotropy in the exchange term 
comes from the orbital ordering with 3 z 1 2 3 4 5 6 7 8 9 10 — r 2 -type or¬ 
bitals pointing in the [111] direction taken as the z axis 


in the calculation 32 . The ratio of J\\/J± was determined 
from a mean-field analysis of the low-temperature mag¬ 
netization curves, with the value of D/J± fixed accord¬ 
ing to the observed cycloidal pitch, a cyc . As long as the 
magnetic states of slow spatial variation such as cycloid 
and skyrmion are concerned, influence from complicated 
crystal structure becomes reduced, which justifies the 
choice to study the phase diagram using a lattice spin 
model on the simple cubic lattice. In the latter case, the 
Dzyaloshinskii-Moriya term is given by 

'Hdm = -D ^2 [y ■ ( m i x rm+x) - X ■ (rrii x m i+ y)\ .(3) 

i 

By comparing the energies of the three magnetic 
phases, we find that only the cycloidal and the ferro¬ 
magnetic states are stable at T=0, while the SkL state is 
not. The transition points at finite T are determined ac¬ 
cording to the calculated temperature-dependence of the 
specific heat and the net magnetization. The specific heat 
has a peak at the transition from the paramagnetic to the 
ordered phases. At the transition from the ferromagnetic 
phase to the cycloidal or the SkL phase, the magneti¬ 
zation shows a sharp drop with increasing temperature, 
while the specific heat shows almost no anomaly. Mag¬ 
netic structures in these phases are identified from both 
the real-space and the reciprocal space mappings of their 
local magnetization configurations. 
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